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We have observed cosmic-ray electrons from 10 GeV to 800 GeV by a long duration 
balloon flight using Polar Patrol Balloon (PPB) in Antarctica. The observation was 
carried out for 13 days at an average altitude of 35 km in January 2004. The detector 
is an imaging calorimeter composed of scintillating-fiber belts and plastic scintilla- 
tors inserted between lead plates with 9 radiation lengths. The performance of the 
detector has been confirmed by the CERN-SPS beam test and also investigated by 
Monte-Carlo simulations. New telemetry system using a commercial satellite of Irid- 
ium, power supply by solar batteries, and automatic level control using CPU have 
successfully been developed and operated during the fiight. From the long duration 
balloon observations, we derived the energy spectrum of cosmic-ray electrons in the 
energy range from 100 GeV to 800 GeV. In addition, for the first time we derived 
the electron arrival directions above 100 GeV, which is consistent with the isotropic 
distribution. 
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1 Introduction 

High-energy cosmic-ray electrons cannot propagate far from the sources, be- 
cause the electrons lose rapidly energy with an energy loss rate of the square 
of energy through the synchrotron and inverse Compton process. Kobayashi 
et al. (2004) suggests that the energy spectrum of cosmic-ray electrons 
might have a structure in the energy region over several lOOGeV due to 
the discrete effect of local sources. This means that we can identify cosmic- 
ray electron sources from the electron spectrum over several lOOGeV to- 
gether with anisotropics in arrival direction. In addition, there is a possibility 
that Weakly Interacting Massive Particles (WIMPs) annihilate directly into 
electron-positron pairs and produce mono-energetic electrons and positrons 
[2]. Although the propagation through the Galaxy would broaden the line 
spectrum, the observed electron and positron spectrum could still have a dis- 
tinctive feature. Since there are no other known production mechanisms that 
would produce an electron and positron peak at energies of 100 GeV - 10 TeV, 
such a distinctive feature clearly indicates the existence of WIMP dark matter 
in the Galactic halo. Thus, the observations of high-energy electrons bring us 
unique information about sources and propagation of cosmic rays, and enable 
us to search for WIMP dark matter. 

Although the cosmic-ray electrons have been observed with many kinds of 
detectors borne on balloons, most observations are limited below ~100 GeV. 
As an exceptional instance, the observation with the emulsion chambers (EGG) 
has achieved to detect electrons beyond ITeV by the large amount of exposures 
of ~ 14 days in total for 30 years [3lH] . The difficulty of observation originates 
from that the electron flux itself is very small and decreases with energy much 
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more rapidly than that of protons because of the electro-magnetic energy loss 
by radiation. The electron flux is estimated to be ~1% of protons around 
lOGeV and less than ~0.2% of protons above 100 GeV from the observed 
energy spectra with a power-law index of —3.0 ~ —3.3 for electrons [5P] and 
—2.7 for protons [5j. Therefore, for the electron observations above 100 GeV, 
we need a long duration observation by a detector with an excellent capability 
of large geometrical factor and powerful background rejection. 

In order to meet the requirements for the electron observation, we have de- 
veloped a highly granulated imaging calorimeter, the Balloon-borne Electron 
Telescope with Scintillating fibers (BETS), that preserves the superior quali- 
ties of both electronic detectors and emulsion chambers. We have successfully 
observed cosmic-ray electrons in the energy range from 10 to 100 GeV at San- 
riku in Japan [6] . The BETS was improved to observe atmospheric gamma rays 
at mountain and balloon altitudes for calibrating the atmospheric neutrino 
flux calculations Furthermore, for observing the electrons over lOOGeV by 
a long duration balloon experiment, we have achieved the development of an 
advanced detector, PPB-BETS, flown by PPB in Antarctica. The PPB has a 
capability to realize a flight for 2 — 4 weeks in one round of the Antarctica at 
an altitude of ~35 km. The PPB project is carried out by the inter-university 
collaboration with Institute of Space and Astronautical Science (ISAS), JAXA 
and National Institute of Polar Research (NIPR) 

In this paper, we describe the instrumentation of PPB-BETS, the calibration 
of the detector with accelerator beams, the observations in Antarctica, and 
the results of our experiments. 



2 Instrumentation 

2.1 Detector 

The PPB-BETS detector consists of 36 scintillating fiber belts, 9 plastic scin- 
tillators, and 14 lead plates with 9 radiation lengths (r.l.) in total. Each fiber 
belt is composed of 280 fibers with a 1mm square cross section for each. A 
schematic cross section of the detector is shown in Fig. [H The area of the 
detector is 28 cmx28 cm and the total thickness is 23cm, including spacers 
inserted between sensitive layers. The scintillating fibers perform the function 
to detect the shower particles developing in lead plates and to fulfill a detailed 
imaging of the shower. Scintillating fiber belts are set in right angle alternately 
to observe the projected shower profile in x and y directions. For the read-out 
of scintillation light in the fibers, we used an image-intensified CCD camera in 
each direction of x and y. The plastic scintillators were adopted for the instru- 
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ment trigger and the energy measurement. Thus PPB-BETS is an instrument 
to observe the details of three dimensional shower development with a timing 
capability. 

The basic structure of PPB-BETS is similar to the BETS [TD]. In order to 
achieve the observation of electrons up to ITeV by using a long duration 
flight, however, the PPB-BETS has been improved in comparison with the 
BETS at several points presented in the following. The thickness of lead is 
increased from 7 to 9 r.l. to observe higher energy electrons, and the number 
of plastic scintillators is increased from 3 to 9 to detect accurately the shower 
development. The image-intensified CCD camera (I.I.-CCD) system for the 
read-out of scintillating fibers has newly been developed by using a better 
quality of CCD to detect ITeV electron showers without saturation of CCD 
signals at higher intensities of photon numbers. 

The total performance of detector was approved by the beam test using ac- 
celerators and by the Monte-Carlo simulations. Monte-Carlo simulations were 
carried out with several codes such as FLUKA2002 pT][T2] . GEANT3 and 
EPICS [H], in which we could not find serious differences. 

2.2 Trigger system 

The trigger rate in actual observation is critically important to know the data 
amount transmitted by telemetry, since we have no guarantee of the recovery 
of instrument. Therefore, we carried out a test flight of the trigger system at 
Sanriku in advance of the flight in Antarctica. The trigger system is composed 
of the plastic scintillators and lead plates in same configuration of the PPB- 
BETS. The flight has confirmed that the trigger system worked well and the 
rate observed was consistent with the expectation by the simulation. We made 
also a study on the detailed performance of detector at CERN-SPS by using 
the electron beams from lOGeV to 200GeV and proton beams from 150 to 
350GeV. The electron events developed purely in electro-magnetic process 
have similar and smooth shower developments. On the other hand, the proton 
events caused by nuclear interactions have the developments with much wider 
fluctuations and relatively smaller energy depositions. The detection efficiency 
of the shower trigger system was examined by adjusting the discrimination 
levels at each scintillator. 

2.3 Imaging system 

For the image-intensified CCD camera, it is required to observe shower images 
from a single track to a shower maximum without saturation. The CCD (THX 
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7887A) has 1024x1024 pixels, and the signal-to- noise ratio is 70dB that means 
6.3 times larger dynamic range of CCD than the BETS of 54dB. Since the 
pixel number increases by 16 times comparing to that used in BETS, the 
dynamic range of CCD for the photon intensity in one fiber extends by two 
order of magnitude. Extrapolating the data of the CERN beam test to higher 
energies, we confirmed that the 1. 1. -CCD system has an enough dynamic range 
for the observation of one single track to ITeV electrons with linearity. The 
one CCD image is read out with a speed of l/30sec by using 12 bit ADC of 
40MHz. Pile-up of events in CCD are taken place only by successive events 
within the timing of 10/isec, and the effect can be negligible small. The total 
data processing on-board take ~ Isec for one event. 

Since the instrument is not scheduled to be recovered, all of the observed data 
should be transfered by the Iridium satellite phone line and the down link to 
the Syowa station. To meet the telemetry rate of 2.4kbps by Iridium, the image 
data can be compressed on-board with a Run-Length method and a Huffman 
method. From the beam test, we obtained that the compressed image data size 
of one event is ~18kByte for electrons at 100 GeV and ~13kByte for protons 
at 100 GeV with the data compression rate of ~40%. 

2.J^ Ballooning technology 

As a ballooning technology of the long duration flight, we have developed 
the telemetry system with the Iridium satellite phone system, power supply 
system by solar batteries, automatic level control system using CPU, and so 
on. Since the detector weight should be saved for loading the heavy ballast, 
we used an un-pressurized vessel with a light shield, designing the instrument 
to meet the vacuum and heat conditions during the long duration flight. We 
examined the validity of the instrument by environmental tests as long as 
more than one week. 

We summarize the basic parameters of the PPB-BETS instrument in Table [TJ 



3 Observations 

The balloon was launched at the Syowa station (69°00' latitude south, 35°35' 
longitude east) in Antarctica at 15:57 on January 4, 2004 (UTC). The level 
flight was started at 18:00 on January 4, and continued till 1:46 on January 
17, 2004 (UTC) at altitudes of 33-37 km (35km on average). The total dura- 
tion of the exposure time is 2.96x10^ hr. The trajectory of the PPB-BETS is 
presented in Fig. [21 The balloon went round the Antarctica at ~65° latitude 
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south from east to west with a speed of ~ 30 — 35 km h~^. Power consump- 
tion of TOW in the instrument was normally supplied by the solar batteries. 
Automatic level control system successfully operated as presented in Fig. [HJ 

The event trigger was executed by two modes, the high-energy (HE) mode 
and low-energy (LE) mode. The LE mode corresponds to electron observa- 
tion over lOGeV, and was assigned for the observation during 10 hours just 
after launching. The data acquired by the LE mode were directly transfered 
to the Syowa station with a telemetry of 64kbps. The HE mode, which cor- 
responds to the electron observation over lOOGeV, was used through all the 
flight. The acquired data were further selected by the software trigger (2nd 
trigger) on-board, and transmitted to our operation room at NIPR in Japan, 
via a receiving station in US, with an Iridium satellite phone line at a rate 
of 2.4kbps. The commands from the operation room to the PPB-BETS were 
also sent by the Iridium line. The monitoring data were delivered to each in- 
stitute via Internet. The data transfer rates were ~ 1x10^ events per hour 
(3Hz) for the LE mode and ~ 70 events per hour (0.02Hz) for the HE mode, 
respectively. The amount of data used for the current analysis is ~ 1/4 of the 
acquired data on-board because of the problem of data transfer system and 
CCD noise. 



4 Data analysis 

4-1 Data sets for analysis 

For the event triggering the HE mode, we adopted two main sets of discrimina- 
tion levels of pulse heights in all these plastic scintillators. One set corresponds 
to the threshold energy of 100 GeV and the other to 150 GeV. The number 
of acquired events is 3.1x10^ for the 100 GeV threshold and 1.6x10^ for the 
150 GeV threshold. There are also 1.0x10'^ events for the other thresholds 
on the HE mode. In this analysis, we used the two data sets of the 100 GeV 
and 150 GeV threshold without use of the other data including the LE mode. 
Figure H] shows the pulse height (converted to the number of Minimum Ion- 
izing Particle, MIP) distribution of the plastic scintillator at the depth of 7 
r.l. for the two thresholds. As shown in Fig. HI the observable minimum en- 
ergy is shifted properly depending on each threshold, and the trigger system 
worked normally. The spectra in high-energy region are compared to a power- 
law spectrum with an index of —2.7. They are consistent with the spectral 
shape of the cosmic-ray protons, since the events triggered on-board are still 
dominated by the protons that are much more abundant than the electrons. 



6 



4-2 Energy measurement 



The plastic scintillators are used for the energy measurements of electrons. 
Since the number of shower particles at the shower maximum is nearly pro- 
portional to the energy of incident electrons, we determined the electron en- 
ergies by using the number of particles at the shower maximum in the transi- 
tion curve with plastic scintillators. The transition curves of electro-magnetic 
shower are well represented by the following formula: 



where S is the pulse heights (converted to MIP) of plastic scintillators at the 
depth of t r.L, N is the normalization factor, 6 is the zenith angle of the shower 
axis, and a and b are adjustable coefficients. In order to derive the number of 
particles at the shower maximum, we fitted this formula to the pulse heights of 
the six scintillators at depths of 3 r.l. to 9 r.l. As shown in Fig. [5], the relation 
of the number of particles at the shower maximum and the energies of electron 
beams at CERN-SPS is almost linear in the energy range of 10 — 200 GeV. 
The energy resolution is 12 % at 100 GeV as shown in Fig. [6|, and consistent 
with the simulations. 

4-3 Angular measurement 

It is important to determine accurately a shower axis, since it is crucial to 
determine the electron energy and the separation of electrons from background 
protons and gamma rays as mentioned later. 

For the acquired events in the fiight, we reconstructed the raw CCD images of 
showers to the fiber positions in detector space by using the positions of each 
fiber on the CCD image. The positions of each fiber were obtained in advance 
of the fiight by observing cosmic-ray muon tracks at ground level. The relative 
displacement of the fiber positions during the fiight was calibrated by using the 
LED-illuminated fibers. Examples of the reconstructed shower image observed 
in the fiight are presented in Fig. [71 

From the reconstructed shower image, we determine the shower axis in the 
following. At first, below 3 r.l. in the detector, we derive centric positions of 
shower lateral spreads in each fiber layer. By fitting these centric positions 
with a linear least-squares method we determine a tentative shower axis, and 
trace back to the next upper layer. Adding a new centric position in the upper 
layer by using only fiber signals near the tentative shower axis (±5 fibers), we 
fit the shower axis once again. This process is repeated up to the top layer of 
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r.L, and the shower axis is determined. Figure [8] shows the angular resolution 
by the beam test, which ranges from 0.40 to 0.57 degree with electron energies 
of 10 to 200 GeV. Although they are considerably better than those of the 
BETS of ~ 1°, they are worse than the simulated angular resolution of ~ 0.15° 
above 100 GeV. It is caused by errors of the fiber positions in detector space. 



4-4 Electron selection 



The electron selection from backgrounds has been carried out by the method 
used in BETS as explained below. At first, proton backgrounds are reduced 
with the on-board trigger system. The discrimination levels are set with each 
scintillator. Second, we selected the events whose shower axes pass through 
the top and bottom of the detector within a zenith angle of 30 degrees. This 
selection reduces the proton background events triggered on-board, which are 
incident from the side of the detector. Third, we introduced the ratio of en- 
ergy deposition within 5mm from the shower axis to the total {RE parame- 
ter) as described in Torii et al. (2001) |6J. A typical event of electron-induced 
shower presents a narrower lateral spread concentrated along the shower axis, 
as shown in Fig. [3 On the contrary, that of proton-induced shower shows a 
wider lateral spread. We characterize this physical property by the RE pa- 
rameter. The electron candidates were selected by the RE distribution under 
the condition that the RE values are greater than 0.75 [T5] . 

As for the separation between electrons and gamma rays, electrons could be 
identified by the presence of hits in the scintillating fibers at r.l. along the 
shower axis. Incident electrons leave signals on the fibers along the shower 
axis. On the other hand, incident gamma rays leave no signals except for the 
back-scattered particles. We judged electrons as the events whose distances 
are less than or equal to 5 mm. Gamma-ray events are rejected by 90 % [15J. 

As a result, the number of the electron candidates above 100 GeV is 84 events 
among the 4.7x10'^ analyzed events in the HE mode. 



4-5 Derivation of electron energy spectrum 



From the electron events, with the above selection and energy determination, 
we derived the cosmic-ray electron spectrum by using the following formula: 

Je{E) = (^g^^^j^^^^cnh — C'2nd)C'atm, (2) 
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where Ne is the number of electron candidates, SQ the effective geometrical 
factor, T the observed live time, AE the energy interval, Cre the correction 
factor of the proton contamination in the RE-cut with energy dependence, 
Ceg the correction factor of the gamma-ray contamination in the gamma- 
ray rejection with energy dependence, Ccnh the correction of enhancement of 
flux due to the energy resolution, C2nd the flux of secondary (atmospheric) 
electrons at the observation level, and C^tm the correction factor of energy 
loss of primary electrons in the overlying atmosphere. 

The effective geometrical factor, SQ, is derived by Monte-Carlo simulations 
under the same condition of the experiment. Figure [3 shows the effective geo- 
metrical factor with electron energies for the 100 GeV threshold and 150 GeV 
threshold in the HE mode. 

The uncertainty of the energy determination has the effect of enhancing the 
absolute flux of electrons, in particular, for a steep power-law spectrum. We 
derived the enhancement factor Cgnh above 100 GeV of 0.98 due to the energy 
resolution presented in Fig. O 

In order to derive the primary electron spectrum, we subtracted the secondary 
electrons produced by the interactions of cosmic-ray nuclei. The energy spec- 
trum of the estimated atmospheric electrons are represented by 

C2ndiE) = l.32xlO-\—§—)-'-'H^-'s-hT-'Gey-') (3) 



at the altitude of 7.4 g cm ^ [31123] . 

The atmospheric correction factor, Catm, is caused by bremsstrahlung energy 
losses of the primary electrons. In the case of a single power-law spectrum with 
an index of —7, the energy-loss correction due to the overlying atmosphere is 
given by a = e^'-'^"^-'*/*^'^^^'', where t is the depth in r.L, and the notation A 
refers to the formula used in electro- magnetic shower theory [3]. The energy- 
loss correction, a, corresponds to the correction factor of flux as follows: Catm = 
a'^~^. In the case of 7 = 3.0 and 7.4 g cm~^ depth, the atmospheric correction 
factor is Catm = 1.37 for a = 1.17. 

From the raw numbers of electron candidates, we evaluated the electron num- 
bers by correcting the proton contamination in the RE cut and the gamma- 
ray contamination in the gamma-ray rejection. Since the power-law index of 
the primary proton spectrum is —2.7, which is much harder than that of the 
electrons, the contamination of protons increases with energy. Atmospheric 
gamma rays are also produced by the interaction of primary protons with at- 
mospheric nuclei. Therefore, the index of gamma rays is the same with that 
of protons. This brings the same energy dependence for the contamination 
of gamma rays with the protons. Hence, the correction factors are estimated 
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to be Cue = 1 — 0.325 x (E'/lOOGeV)'''"^'^ for the proton contamination and 
Cgg = 1 — 0.176 X (i?/100GeV)'*'~^''^ for the gamma-ray contamination from 
the Monte-Carlo simulations. 

4-6 Check on event identification 

In order to check on the event identification, we derived the zenith angle 
distribution of electron and gamma-ray events. For cosmic-ray electrons, the 
zenith angle distribution will be isotropic, and thus for a plane detector the 
observed angular distribution will be proportional to cos9d{cos9), neglecting 
the effect of the finite detector seize. The distribution for atmospheric gamma 
rays should be constant with cos6', since the gamma rays are produced in the 
overlying residual atmosphere and thus have a 1/ cos^ enhancement relative to 
the isotropic primary hadron fiux. Zenith angle distributions for the electrons 
and gamma rays are presented in Fig. [T0|, and in good agreement with the 
expected forms with the finite detector size. 

We also derived the vertical spectrum of atmospheric gamma rays at 7.4 g cm~^ 
originated from hadronic interactions, subtracting the bremsstrahlung gamma 
rays produced by primary cosmic-ray electrons. Figure [11] shows the observed 
spectrum of atmospheric gamma rays at 7.4 g cm~^, which is well represented 
by 

= (3.19±0.59)xl0-^(^^^|-^)-2-7^±°-2^(m-Visr-iGeV'i). (4) 

The spectral shape is in good agreement with the ECC and the calculated 
results from hadronic interactions of primary cosmic rays with the overlying 
atmosphere [23] • At the normalized depth of 7.4 g cm~^, the gamma-ray fiux 
agrees with the BETS in the energy region of 10 - 20 GeV [7] , and is consistent 
with that of ECC within 95% confidence level of statistical errors. 



5 Results and discussions 

5. 1 Energy spectrum of electrons 

Figure [12] presents the electron energy spectrum multiplied by the cube of 
energy, in comparison with other electron observations [ r6|17f 18|[T9][6][20ll^][2T] . 
Table [2] presents the fiux of electrons at the top of atmosphere in each energy 
interval. Our spectrum agrees well with the extrapolated spectrum of the 
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BETS [6] at 100 GeV. The combined spectrum of PPB-BETS and BETS can 
be represented by a single power-law function of 

J, = (1.82±0.13)xl0-^(^^^|^)-3-™5(m-2s-isr-iGeV-i) (5) 



in the energy range of 10 GeV to 800 GeV. The reduced chi-square value is 
1.605 for 11 degrees of freedom. A power-law spectrum is acceptable at the 
95 % confidence level of statistical errors. 

The energy spectrum exceeding 100 GeV is crucial to detect the nearby SNRs 
as discussed by Kobayashi et al. (2004) [1], and electron-positron pairs from 
Kaluza-Klein dark matter annihilations [2]. Although the data statistics of 
our results are insufficient to discuss the details of the contribution of nearby 
SNRs and/or WIMP dark matter, our energy spectrum may indicate a sign of 
a structure in the several 100 GeV region, as shown in Fig. [121 Similar structure 
in the electron energy spectrum is reported by the ATIC-2 observations |21j . 
Therefore, although the energy spectrum with PPB-BETS cannot reject a 
single power-law function, both the observations with PPB-BETS and ATIC- 
2 may indicate a significant spectral structure. For the future observations, we 
are planning to increase greatly our statistics by experiments such as CALET 
with the geometrical factor of nearly 1 m^sr and the observation time of 3 years 
on the International Space Station [22]. It is expected that CALET can detect 
these distinctive signatures from the nearby SNRs and WIMP annihilations 



5.2 Arrival directions of electrons 



Incident directions of electrons on the PPB-BETS detector are determined by 
using the shower axis with an accuracy of around 0.5 degree. Since the atti- 
tudes of PPB-BETS instrument are determined with sun aspect sensors and 
geomagnetic aspect sensors, we can determine arrival directions of cosmic-ray 
electrons with an accuracy of several degree. For the determination of attitudes 
with geomagnetic aspect sensors, we referred to the International Geomagnetic 
Reference Field (IGRF) model j24|. From these capabilities of PPB-BETS, for 
the first time we derived the arrival directions of high-energy electrons. Fig- 
ure [13] presents a ratio of the observed distribution above ~100 GeV to the 
isotropic distribution along the Galactic longitude. As shown in Fig. [131 the 
arrival directions of electrons are consistent with the isotropic distribution. 

Because the rate of energy loss due to the synchrotron and inverse Compton 
processes is much higher for electrons than for nuclei, the degree of anisotropy 
of high-energy cosmic ray electrons is expected to be higher than that of the 
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nuclear component [251I26] . The magnitude of anisotropy is expected to be 
only around 1 % in the 100 GeV - 1 TeV in the direction of Vela at (£, b) = 
(— 96°.l, —3°. 3) by using the calculated results by Kobayoashi et al. (2004) 
P12S] . The first result of the electron arrival directions by PPB-BETS shows 
no significant anisotropy within statistical errors. For the future observations, 
the detection of anisotropy toward the sources will enable us to identify the 
cosmic-ray electron sources, together with the unique signatures of the energy 
spectrum described above. 
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Table 1 

Basic parameters of the PPB-BETS 



Tn ti'I'yn tyi cini" \KTC•^ crh'l' 




mr'niriTncr nn— t^tpqcitti ^(^r\ o"An H a1 5i 


1 Tlnf a 1 AST'pi o"Vi f 




inplnrlincr Hftllj^cjf 1 


AurPT p An cii 1 TY~i T\i' 1 An 


7f) W 

1 \J V V 


t^mrnliprl '^Alfir Hat I'PTipt? 


Level altitude 


^SS km 


by automatic level control 


Data transfer rate 


2.4 kbps 


by the Iridium phone line 




64 kbps 


by the down link to the station 


Energy range 


10 GeV-1 TeV 


by two trigger modes 


Effective geometrical factor 


~300 cm^sr 


> 100 GeV by simulation 


Energy resolution 


10-20 % 


in 10 GeV-1 TeV by plastic scintillators 


Angular resolution 


0.4-0.6 deg 


by scintillating fibers 



Table 2 

Electron flux at the top of atmosphere 



E 


E 


Intensity" 


(GeV) 


(GeV) 


(m-^s-^sr-^GeV-i) 


100.0-158.5 


122.6 


(8.1±1.6)xl0-5 


158.5-251.2 


194.4 


(3.4±0.7)xl0-5 


251.2-398.1 


308.0 


(9.0±2.4)xl0-^ 


398.1-631.0 


488.2 


(2.8±0.8)xl0-^ 


631.0-1000.0 


773.7 


(2.0±0.9)xl0-^ 



" The quoted errors are the statistical error. 
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Fig. 2. The trajectory of the PPB-BETS in Antarctica. 
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Fig. 4. Pulse height distributions of the plastic scintillator at the depth of 7 r.l. 
for the 100 GeV threshold (solid circles) and 150 GeV threshold (open circles). The 
solid line shows a power-law spectrum with an index of —2.7. 
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Fig. 5. Relation of the electron energy and the number of shower particles at the 
shower maximum by the CERN-SPS beam tests with simulations. The incident 
direction is perpendicular to the detector surface. 
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Fig. 6. Energy dependence of the energy resolution for the CERN-SPS electron 
beams with simulations. The incident direction is perpendicular to the detector 
surface. 
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Fig. 7. Examples of the reconstructed images in detector space. Typical candidates 
of an electron, gamma ray, and proton are shown in x and y directions. 
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Fig. 8. Energy dependence of the angular resolution for the perpendicular electron 
beams at CERN-SPS. 
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Fig. 9. The effective geometrical factor SO, with the Monte-Carlo simulations under 
the same condition of the experiment. 
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Fig. 10. Zenith angle distributions for events classified as electrons and gamma rays 
with the expected distributions. 
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Fig. 11. Energy spectrum of atmospheric gamma rays observed with PPB-BETS, 
compared to the BETS [7J and ECC [23j. The gamma-ray fluxes are normahzed to 
7.4 g cm~^ equivalent altitude. The dash line shows the best flt power-law function 
of PPB-BETS data with an index of -2.77±0.21. 
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Fig. 12. Electron energy spectrum observed with PPB-BETS (solid circles) in com- 
parison with the energy spectra of BETS (solid squares) ant the other observations. 
The dash line shows the best fit power-law function of the combined spectrum of 
PPB-BETS and BETS with an index of -3.05±0.05. 
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Fig. 13. A ratio of the electron flux observed with PPB-BETS to the isotropic 
distribution along the Galactic longitude. 
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